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ABSTRACT 

We argue that current observations of M — a relations for galaxies can be used to constrain theories 
of super-massive black holes (SMBH) feeding. In particular, assuming that SMBH mass is limited 
only by the feedback on the gas that feeds it, we show that SMBHs fed via a planar galaxy scale 
gas flow, such as a disc or a bar, should be much more massive than their counterparts fed by 
quasi-spherical inflows. This follows from the relative inefficiency of AGN feedback on a flattened 
inflow. We find that even under the most optimistic conditions for SMBH feedback on flattened 
inflows, the mass at which the SMBH expels the gas disc and terminates its own growth is a factor 
of several higher than the one established for quasi-spherical inflows. Any beaming of feedback away 
from the disc and any disc self-shadowing strengthens this result further. Contrary to this theoretical 
expectation, recent observations have shown that SMBH in pseudobulge galaxies (which are associated 
with barred galaxies) are typically under- rather than over-massive when compared with their classical 
bulge counterparts at a fixed value of a. We conclude from this that SMBHs are not fed by large (100 
pc to many kpc) scale gas discs or bars, most likely because such planar flows are turned into stars 
too efficiently to allow any SMBH growth. Based on this and other related observational evidence, 
we argue that most SMBHs grow by chaotic accretion of gas clouds with a small and nearly randomly 
distributed direction of angular momentum. 
Subject headings: galaxies: formation - galaxies: active - accretion: accretion discs 



1. INTRODUCTION 

It is now well established that super-massive black 
holes (SMBHs) reside in the nuclei of many galax- 
ies. The masses of these SMBHs correlate with a 
range of properties of their host sphero ids, includ- 
ing luminosity (e. g. iMagorrian et all Il998|) and conse- 
quen tly mass (e.g. iMarconi fe Huntl l2003: Hii ring fe Rixl 
20041): the concentr ation of stellar bulge light (e.g. 
Graham et al.1 120011 ): the deficit o f stellar bulge 



light i n ellipticals (e.g. iGrahaml l2004t iFerrarese et al.l 



120061: IKormendv fe Benderl I2009D : and the gravita- 
tiona l binding energy of t he bu lg e (e.g. IFeoli fe Meld 
2001 Ide Francesco et all 12001 IFeoli fc Meld 120071: 



Aller fc Richstonell2007D . 



The correlation between SMBH mass (Mbh) and the 
velocity dispersion (a) of its host spheroid, referred 
to below as a "Mbh-c relation", has been studied by 
many authors. A power-law fit to the data, Mbh oc 
<7 V , y ielded values of p in th e range of p ~ 4 — 
5 Ccf. IFerrarese fc Merrittl 120001 iGebhardt et all |2000l 
iTremaine et all 120021 iGultekin et all 120091 ). and is now 
widely interpreted as the imprint of the self-regulated 
growth of the SMBH. A more recent extended sample 
of galaxies sho ws that the power la w index p may be as 
high as p w 6 (| Graham et al.ll2011[) . 

Furthermore, recent work has highlighted that the 
Mbh _ c relation one measures will depe nd on the nature 
of the galaxy sample. For example, iMcConnell et all 
(J2011I ) presents two extremely massive SMBH, Mbh ~ 
1O 1O M0, in the centres of two giant elliptical galaxies. 
Combined with other previous SMBH mass measure- 
ments, these authors obtain p rj 4.5 for elliptical and 
SO galaxies. A similar but offset to lower SMBH masses 



result is found for spiral galaxies in the sample. Alter- 
natively, when data for all galaxies are combined in one 
Mbh _ cr relation, p rs 5 is found. 

Another important conclusion is gradually emerging 
from the observations: SMBHs in galaxies that show flat- 
tened morphology, such as discs, pseudobulges or bars 
appear to be less massive at same a than their coun- 
terparts in elliptical galaxies. In particular, iHul (|2008t ) 
showed that pseudo bulge galaxies t end to have under- 
massive SMBHs, and I Graham! (|200 8) showed that barred 
galaxies tend to have SMBHs that are a factor of ~ 3 to 
ten smaller than el liptical galaxies at the same velocity 
dispersion. While IKormendv et al.l (J2011I ). in variance 
with other authors, find that SMBHs in pseudobulges 
and galaxy discs do not even follow an Mbh _ cr relation, 
their results are consistent with that of the other authors 
in that most of their ps eudobulge sys t ems li e below the 
Mbh-f for ellipticals. iMathur et al.l (|2011| ) firmed up 
these conclusions with a sample of narrow line Seyfert 1 
nuclei (all classifiable as pseudobulge systems). A useful 
compilation of the recent data, summarisi ng these points, 
is give n in the right panel of Fig. 1 in iShankar et all 
(12011 . 

The purpose of our paper is to show that there is a uni- 
fied theoretical picture within which this observational 
trend can be understood, and that the observed Mbh" 
a relations place very interesting constraints on the still 
unknown feeding mode of SMBHs. 

A well known difficulty in fuelling active galactic nu- 
clei (AGN) is that the typical angular momentum of gas 
in the galactic bulge is very large compared wit h that 
of the last stable orbit around a black hole (e.g., iKrolikl 
119991: lCombesll200ll: Uogeell2004D . One expects that due 
to circularisation shocks, gas will end up in a planar fea- 



ture, such as a disc or a ring (cf. simulation S30 shown in 
Figur es 1 and the left panel of Figure 13 in lHobbs et al.l 
1 2 llf ) . It has been argued that large scale angular mo- 
mentum transfer occurs rapidly enough through the ac- 
tion of stellar and gaseous bar s (e.e., [S hlosman et a . 
I1990T ) or via spiral density waves (jThompson et al.1l2005f l 
to overcome this difficulty. These galaxy-wide gas discs 
(with radial scales of many kpc) could then extend all the 
way to the inner galaxy and even feed the SMBH. This 
"planar mode" of feeding is schematically illustrated in 
the left panel of Figure 1. 

However, gas discs tend to become unstable to 
gravitational fragmentation and resulting star forma- 
tion (S F) beyond a self-gravity radius of ~ 0.1 
dc (e.g.. iPaczvhskil [19781 iKolvkhalov fe Sunyaevl J l980l: 
iShlosman fe Begelmanl 119891 : iCollin fe Zahnl 11991 ). In 
the simplest model this would result in complete 
consumption of gas i n the disc by star formation 
(|Navakshin et al.l |2007|) and so there would be no way 
to fuel the SMBH. However, stellar feedback could heat 
up the large scale disc to self-regulate the star formation 
rate in the disc to a level such that a sufficient amount 
of gas can reach the S MBH, but there is no clear consen- 
sus o n this point (e.g.. lGoodmanll2~003l : lThompson et al.l 
f2005h . 

An alternative solution to the diffic ulty of fuelling 
AGN is provided by chaotic feeding ([King fe Pringld 
120071 ) or ballistic accretion of cold streams when the 
gas is cold enough to form dense clouds or filaments 
(jNavakshin fc Kind [20071 : IHobbs et al.l l20ll . Large 
scale cosmological simulations show that cold streams 
may be important for supplying rel atively low angu- 
lar moment um gas to galaxies(e.g., iKeres et"aT1 120091 : 
iKimm et afll201l[ iDubois et al~ll2011[ ) , while higher res- 
olution simulations that resolve the regions nearer to 
SMBH suggest that similar proce sses may operate on 
smaller scales ([Levine et al.l 12010ft. In part icular, the 
numerical experiments of IHobbs et all ([2011D show that 
turbulence (driven by stellar feedback in the bulge) in a 
quasi-spherical distribution of gas leads to formation of 
convergent flows that create high density filaments. The 
latter can travel almost ballistically through the rest of 
the bulge, with some filaments arriving in the inner par- 
sees of the galaxy. Note that most of the gas still ends 
up in a galaxy-scale disc with a large angular momen- 
tum; however, the SMBH receives most of its fuel not 
from that disc but from the gas filaments arriving at its 
vicinity directly. 

This "chaotic cold streams" feeding mode is schemati- 
cally shown in the right panel of Figure 1. In this class of 
model the orientation of the innermost sub-pc scale disc 
makes a random walk because of the arrival of material 
with a fluctuating angular momentum. The cancellation 
of the angular momentum in shocks transfers gas into 
the SMBH more rapidly than viscous torques can, and 
so these discs should be more resilient to star formation. 
Furthermore, because the orientation of the inner disc 
is constantly changing, feedback from the AGN will be 
quasi-isotropic. 



The main point of our paper can be summarised 
as follows. Large scale gas discs are better able to 
withstand AGN feedback than quasi-spherical inflows 
(|Navakshin fe P ower 20 13). This is so even under the 



most optimistic assumptions about AGN feedback - that 
is, feedback acting directly on the disc without being 
shadowed by the inner disc, and no beaming away from 
the disc/galaxy plane. We show that, under these most 
favourable of conditions for which AGN feedback is max- 
imal, there is a critical SMBH mass at which its feedback 
expels the disc that feeds it. We call this a "disc Mbh~ 
cr relation" because it arises only for SMBHs accreting 
from large scale galaxy discs. This mass is larger than 
the canonical Mbh-o" relation, and any dilution of AGN 
feedback would increase it even further. Theoretically we 
see that if SMBHs were fed by planar accretion then we 
could expect SMBHs in pscudobulgc systems to lie above 
the Mbh-cr relation for classical bulge s. Observationally 
there is no strong evidence for this (e.g.[Hu 2008; G raham! 
120081 : iKormendv et al.l 120 111 ), and so we conclude that 
most SMBHs do not grow via the planar accretion mode 
but are much more likely t o be fed by quasi-sph erical 
inflows as in the m odels of King fc Pringld ([2007D and 
IHobbs et all (I20T1 . 

To avoid misunderstanding, we note that galaxy-scale 
gas discs and gas bars are obviously important for galaxy 
evolution as a whole, e.g., as a birthplace for a significant 
fraction of all stars, and also by effecting a radial redis- 
tribution of gas and stars within the galaxy, but we see 
no clear observational evidence that these features are 
able to channel their gas all the way to the SMBHs. 



2. THE GAS WEIGHT ARGUMENT FOR A SPHERICAL 
AGN FEEDBACK 

We briefly review existing ideas of AGN feedback act- 
ing on quasi-spherical distributions of gas before mov- 
ing o n to the problem of AGN feedback acting on a gas 
disc. iSilk fe Reesl (J1998I ) assumed feedback in the form 
of an energy-conserving ou tflow and obta i ned a s caling 
Mbh ex cr 5 Subsequently iFabianl (|1999f ). iKingj ([20031 
I2005D and iMurrav et all (|2005l) solved the momentum 
equations for momentum-conserving AGN feedback act- 
ing on a spherical shell of gas in an isothermal bulge po- 
tential (see below) and showed that it leads to a M^-a 
relation 
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where we follow King's analysis; here n is the opacity, 
which is assumed to be dominated by the electron scat- 
tering, and f g is the baryon fraction, which is assumed 
to be equal to th e initial cosmological value of f g ~ 0.16 
fcf. lSperge1l2007f ). We note that all of the aforementioned 
theoretical approaches use d a singular isothermal sp here 
potential (e.g., §4. 3. 3b in iBinnev fc Tremaind [20081 ) for 
simplicity. For such a potential, the one dimensional 
velocity dispersion is a constant independent of radius, 
a = (GAf tot ai/2 J R) 1/2 , where M tota i(R) is the total en- 
closed mass including dark matter inside radius R (the 
distance from the centre of the galaxy) . The enclosed gas 
mass, M(R) = f g M tota \{R) = 2f g a 2 R/G is proportional 
to R. The gas density at radius R for such a potential is 
p (R) = f n a 2 /(2-KG R 2 ) for reference. 
iKingj (|2003l . I2005D assumed that the SMBH luminos- 



ity is limited by the Eddington value LEdd, and the 
momentum outflow rate produced by radiation pressure 



(|King fc Poundsll2003D is of order 
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We can recover the result of IKind (J2003L 120051 ) using 
a simpler "weight argument" , which requires that the 
momentum output produced by the SMBH (IIsmbh) 
should just balance the weight of the gas in the bulge, 
W{R) = GM(R){M totli i(R)]/R 2 , which turns out to be 
independent of radius: 
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(3) 



To order of magnitude, equation |3] holds for any poten- 
tial if estimated at a radius which encloses most of the 
gas potentially available for SMBH fue l ing, e .g., such as 
the virial radius for the iNavarro et al.l (|1997|) potential. 
Balancing the outward force of the outflow (equation [2]) 
with the weight of the gas in the bulge (equation [3]) then 
leads naturally to equation Q] 

The momentum feedback model is attractive in its 
physical simplicity. The model is underpinned by obser- 
vations of fast v ~ 0.1 — 0.3c outflows from bright AGN 
and quasars (e.g., iPounds et alJ 2003a. b). Studying the 
blue shifted ionised absorption features in the 7 to 10 
keV ene rgy range in a l a rge sa mple of AGN with XMM- 
Newton, iTombesi et al.l (|2010f ) concluded that such fast 
AGN outflows appear in ~ 40 % of the sample. As 
these authors emphasise, their uniform sample overcomes 
"publication bias" and shows that fast AGN outflows 
are widespread. It also requires the outflows to be 
wide-angled (rather than jet-like) , which is beneficial in 
spreading the influence of AGN feedback as broadly over 
the galaxy bulge as possible. 

Equation [T] contains no free parameters, but never- 
theless it is very close to the observed M^-a rela- 
tion. Similar reasoning helps to explain both the ob- 
served Mhh-cr relation for Nuclear star Clusters (sec 
iMcLaughlin et all 120061) and the observation that NCs 
are preferentially foun d in low mass (low a) galaxies 
(|Navakshin et al.l I2009T). The model has been used re- 
centlv by IZubovas et al.l (J2011D to explain the two ~10 
kpc scale symmetric lobes apparently filled with cosmic 
rays detected by FE RMI-LAT telescope in the Milky 
Way (|Su et al.1 120101 ) . To explain the particular geom- 
etry ofjhe bubbles, the only adjustment to the b asic 
IKind (120051 ) model required bv lZubovas et al.1 (|2Qlll) has 
been an addition of a dense disc of molecular gas, known 
as the Central Molec ular Zone, found in the central ~ 
200 pc of our Galaxy (jMorris fc Serabvnlll996l ). 

While it is clear that the exact Mbh — cr relation de- 
rived analytically is affected by the choice of the poten- 
tial, we note that the main conclusion of our paper - the 
fact that it is harder to expel a disc feeding the SMBH 
rather than a quasi-spherical shell - remains unchanged 
as it is based on a simple geometrical argument (see be- 
low). Furthermore, gravitational potentials of elliptical 
and lenticular galaxies appear to be relatively well ap- 
proxi mated by a const a nt a value inside their effective 
radii. iCappellari et al.1 (|2006l ) find that the ID velocity 
dispersion of stars within the galaxies in their sample 
depends on radius R (within the galaxy) as a very weak 
power-law: a(R) oc i? -0066 . This implies that a drops 



on average by the factor of ~ 1.16 when R changes by an 
order of ma gnitude, and a A varies by th e facto r of rj 1.8. 

Indeed, IMcQuillin fc McLaughlml (|2012f ) consider 
SMBH momentum feedback for potenti als other than 
isoth er mal, e.g., t he N FW p otential (INavarro et al. 
I1997D iHernquistl (|1990f ) and iDehnen fc McLaughlin 
(|2005l) potentials. In all cases momentum driven out- 
flow produce relations of the form Mbh oc V^, where V^n 
is the maximum circular velocity in the potential. The 
normalisation of the relation is slightly smaller than for 
the simpler isothermal potential we use here, which im- 
plies that escape from the isothermal potential is slightly 
more demanding t han for all the others. To quote a nu- 
merical factor here, IMcQuillin fc McLaughlin] (J2012D find 
a change in the critical SMBH mass (only) by a factor of 
order unity only for the NFW potential, as an example, 
if one defines a at the peak of the rotation curve of the 
potential as a — V m /\/2. 

No analytical treatment for the complex problem at 
hand can be expected to be accurate within better than 
a factor of ~ 2, and therefore we feel that the singular 
isothermal potential is an appropriate approximation to 
use here. 

3. THE THEORETICAL DISC M BH - a RELATION 

As shown bv lNavakshin fc Power! (|2010l ). the above ar- 
gument does not quite work for a non-spherical geometry. 
To see this, consider the simplest example - an axially 
symmetric rotating galaxy. The gas settles in the plane 
of symmetry of the galaxy and forms a rotationally sup- 
ported disc. Let H be the local (i.e., at a given R) verti- 
cal disc scale-height. Let us consider a disc at radius R. 
The maximum momentum flux from the SMBH outflow 
striking the disc is 



H LEdd 

~R c 



(4) 



This is an absolute maximum given by the total momen- 
tum outflow rate (equation 2) times the fraction of the 
solid angle the disc subtends as seen from the SMBH lo- 
cation. In reality some of the momentum outflow from 
the AGN could be intercepted (shadowed) by the disc 
at smaller radii (cf. the left panel of Figure 1), and 
also the outflow may be beamed along the direction per- 
pendicular to the disc. These individual effects actually 
strengthen the conclusion we draw because they act to 
lessen the effect of feedback on the disc and therefore 
should increase the mass of the SMBH. 

We first assume that the disc mass is dominated by 
gas, and we correct for the presence of stars within the 
disc later on. The effective weight of the gas disc is 
given by W e g ~ GM tot ai(i?)M disc /i? 2 , where M disc (i?) = 
ttR 2 T,(R) is the disc mass, and E(iJ) = 2H(R)p c (R) is 
the gas disc surface density, with p c (R) being the mid- 
plane disc density. We note here that a disc in circu- 
lar rotation actually has a zero weight because grav- 
ity is exactly balanced by the centrifugal force. How- 
ever, within a factor of a few the effective weight defined 
above still applicsQ if we wish to expel the disc to infin- 
ity (where the centrifugal force becomes negligible if the 

1 This can be seen by observing that escape velocity from dis- 
tance R to infinity for a circularly rotating disc is smaller by a factor 
of y/2 than that from a static point within the same potential 



angular momentum of the material is conserved). There- 
fore, by requiring that momentum flux balances weight 
(liaise ~ Weff) we limit the SMBH mass to 

k_ M tota \{R)M disc no 2 M disc 
4tt 



M hh = — 



RH 



2nG H 



(5) 



A casual look at this equation would seem to suggest 
that, because H and Mdi SC may vary with radius R and 
from system to system, one should not expect any robust 
correlation between Mbh and a 4 to arise. 

What is missing in the above analysis is the possi- 
bility of star formation in the galaxy disc (this was 
also not included in th e numerical experiments of 
iNayakshin fe Power] 120101 ). Massive cool gas discs be- 
come gravitationally unstable when t he Toomre 
parameter approaches unity from above ()Toomrelll964 



Q 



Kgc s 
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(6) 



where kj! = is the cpicyclic frequency. Noting that kq = 
\/2fl for the singular isothermal sphere potential, Vt 2 = 
GMtotai/R 3 , and using the hydrostatic balance condition 
H = c s f2 _1 , we find that 



Q = 



v/2ffMt ota i(fl) 
RM disc 



(7) 



Note that Mdi sc , defined just below equation HJ is func- 
tion of R in the above equation. 

It is widely believed that self-gravitating gas discs 
self-regulate their star formation rates to maintain 
marginally stable states. As a disc becomes self- 
gravitating, stars form on a local dynamical time scale 
(i.e., extremely rapidly) and the energy and momentum 
released by the stars heat the disc up; this increases Q 
and the star formation rate is reduced to an equilibrium 
rate that maintains quasi-equilibrium. Generally, the 
disc first becomes unstable when Q sa 1.5, and there- 
fore we arrive at the constraint 



H 



Mtotal(R) 

R 



(8) 



Using this now in equation[5l we arrive at the "disc Mbh _ 
a relation" , 

Mbh = ^<x 4 , (9) 

where we introduce /^ < 1 as the gas mass fraction in 
the disc to account for the fact that a realistic galaxy 
disc contains gas and stars as well. According to this 
definition, the fraction 1 — fd of the disc mass is in the 
recently formed stars or very dense massive molecular 
clumps on the way to collapsing and forming stars. 

We see that the disk Mbh~c relation is of exactly same 
form as the "spherical" Mp ^-a relation (equation [1]) de- 
rived bv Kind (|2003U2005[) . However, we note that the 
coefficient in front of the latter is smaller than that in 
the former. The difference is difficult to quantify ex- 
actly, but would seem to be at least a factor of a few to 
ten. Physically, it is due to (1) the factor f g = 0.16 in 
front of the spherical relation and fd for the disk relation. 
It is unlikely that fd is significantly smaller than unity 
during the early gas-rich epoch of galaxy formation. Fur- 
thermore, in principle a similar in meaning coefficient - 



taking into account that part of the spherical shell is al- 
ready converted into stars - may be introduced in the 
spherical M^-cr relation, in which case the difference 
between the two relations would be exactly equal to f g . 
(2) When deriving the disc Mbh — & relation, we took the 
most optimistic assumptions about the AGN feedback on 
the disc. As discussed below equation @] realistically, a 
part of the feedback can be shadowed by the inner disk 
(the innermost part of which is very hard to affect by 
the feedback, see below). Beaming away from the disc 
plane may decrease the momentum feedback on the disc 
further yet. 

Another important point is that although our argu- 
ments were presented for a disc geometry, they are also 
applicable to a bar-mediated gas in flows, such as those 
discussed by iShlosman et al.l (J1990I ) . These authors ar- 
gued that "bars within bars" may effectively channel the 
gas into the central regions of galaxies to fuel AGN. How- 
ever, bars present even less of an effective solid angle than 
discs for AGN feedback to work on: barred gas inflows 
are squashed geometrically not only to the galaxy mid- 
plane, as disc inflows, but also in the azimuthal direc- 
tion within the galaxy plane. Therefore the momentum 
striking a bar would be even less than that for the disc 
(equation 3|) , and thus we would expect even more mas- 
sive SMBHs if they were limited by their feedback on the 
gas bars only. 

The derivation in this section is applicable to regions 
well beyond the black hole's influence radius, Rh = 
GMbh/2c 2 , which is typically of the order of a few to 
tens of pc. In the Appendix we show that a gas disc at 
R <C Rh is actually able to withstand the momentum 
feedback of the AGN. This is why the sketch in Figure 1 
shows small scale discs unaffected by the feedback. 

4. WHY DO OBSERVED SMBHS NOT FOLLOW THE DISC 
Mbh-o" RELATION? 

iKormendv fc Gebhardtl (f200ll ) found that SMBHs do 
not co rrelate with circular velocities of galaxy discs. iHul 
(|2008f ) showed that SMBHs in pscudobulge galaxies are 
underweight by a factor of about ten with respect to 
their cousins in classical bulge systems. iGrahaml (|2008l ) 
showed that barred galaxies tend to have SMBHs that 
are a factor of ~ 3 to ten smalle r than elliptical galaxies 
at the same velocity dispersion. IKormendv et al.l (|2011[ ) 
found that observed SMBHs do not correlate well with 
velocity dispersions of either galaxy discs or pseudob- 
ulges. 

In the sample of IKormendv et al.l ()201lD . only two 
SMBHs in pseudobulge systems appear to lie above the 
classical Mbh-c relation for bulges, with the rest being 
som ewhat belo w the relation, in a qualitative agreement 
with E3 ([111). 

If these observations are confirmed in the future, then 
this is clearly very strongly at odds with the theory de- 
veloped above: because discs are difficult to affect by 
AGN feedback sim ply due to geometrical argu ments (see 
also simulations bv INayakshin fc Pow er 2010), a SMBH 
fed by a disc should be able to grow more massive by at 
least a factor of a few (compare equation O to equation 
[T]) than a SMBH embedded in a bulge. In fact, theo- 
retically one may expect disc SMBHs to be even more 
massive: as noted below equation [4j this is the absolute 
minimum to which the SMBH can grow if fed through 



Planar AGN feeding mode 



Chaotic streams AGN feeding mode 



Galaxy disk 
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Fig. 1. — Schematic illustration of the two modes of AGN feeding. Left - "Planar mode", where gas is delivered from the Galaxy- wide 
(many kpc) disc by a planar inflow through a disc, for example by spiral density waves or "bars within bars". The inner pc-scale disc (light 
blue) is always oriented the same way as the galaxy disc. Right - the "chaotic mode", where the gas is delivered by cold dense filaments 
or blobs on nearly radial orbits. The small scale AGN disc fluctuates in its spin orientation, so that the AGN feedback is quasi-spherical 
when averaged over long (~ million yrs) time scales. The outer Galaxy disc is disconnected from the innermost disc by star formation that 
turns gas into stars before a significant inflow could happen. 



the disc because in general AGN feedback impacting a 
given disc radius can be diluted at smaller radii or be 
beamed away from the disc (neither of these effects is 
included in equation [9]). 

The simplest interpretation of this result is to accept 
that star formation in large scale discs is too efficient a 
process and so significant gas flow into the sub-pc vicin- 
ity of SMBH is suppressed, leading t o a quenching of 
SMBH fueling (see the arguments by [Goodman! [20031 ; 
iSirko fc Goodmanl 12001 iNavakshin et all 120071) . There- 
fore, we conclude that these discs do not extend to the 
SMBH sphere of influence (between 1 and 10 pc, depend- 
ing on SMBH mass), as suggested in the right panel of 
F igure 1. 

iThompson et al.l (|2005| ) proposed a star-burst disc 
mo del which a ppear s to overcome the difficulties noted 
by IGoodmanl (|2003l ). Specifically, they argue that the 
star formation rate in the disc is limited by the action 
of energy and momentum feedback from massive stars 
and supernovae, allowing some fuel to trickle down all 
the way to the nucleu s. However, one problem with 
IThompson et a.1.1 (J2005I ) argument applied to the inner 
galaxy is a mismatch of time scales. Gravitational col- 
lapse of the dis c is expected to take place on the dynam- 
ical time scale (|Gammidl2001l ). 



^dyn 



^3/2 



G l ' 2 M hh 112 



10 3 Rl( 2 M s 1/2 years 



(10) 



where M§ is the SMBH mass i n units of 10 8 M fT ), and 
R pc is radius in units of parsec. IThompson et al.1 (J2005I ) 
show that the disc cooling time is much shorter than the 
dynamical time in their model. On the other hand, mas- 
sive stars dominating the feedback release energy on a 
time scale of a few million years. Therefore, for radii 
where tdyn is much shorter than the lifetime of the mas- 
sive stars, the feedback might be released too slowly, i.e., 
only when the disc would already have collapsed gravi- 
tationally. 

In c o nclusion, we bel ie ve that the observation s of iHul 
(120081) : IGrahaml (120081) : IKormendv et~aT1 (I20TI imply 
that SMBHs do not grow by a large scale (e.g. tens of 
pc to it kpc) discs. 



5. CHAOTIC/BALLISTIC FEEDING OF SMBHS 

Surveys show little evidence for correlation between 
the orientation of jets from AGN and the large-scale 
structure of their host galaxies (jSchmitt et al.1 l2002t 
iVerdoes Kleiin fc de Ze cuw 20051 e.g). This shows that 
the inner £ pc scale gas discs feeding AGN are also 
uncorrelated with the large scale discs of host galaxies. 
In our own Galaxy, t he observed young st ellar discs in 
the Galactic Centre (jPaumard et al.1 120061 ) are inclined 
at very large angles to the Galactic plane. This indicates 
that the disc stars originated from a deposition of one or 
two gas clouds with angular momentum directions very 
different from that of the Galac tic disc (jBonnell fc Ricel 
l2lMlHobbs fc Navakshml[2009l) . 

An attractive physical picture for AGN feeding that 
explains these observations and the results of our paper 
is that SMBHs in general are fuelled by stochastic de- 
position of gas clouds with randomly oriented angular 
momentum, as sketched in t he right panel of Fi gure 1. 

Numerical simulations by iHobbs et al.l ()201lD demon- 
strated that this situation may be realised in an ini- 
tially coherently rotating gas sphere if strong turbulence 
(driven by supernova explosions due to star formation 
in the bulge, for example) is present. Convergent tur- 
bulent gas flows lead to formation of high density re- 
gions travelling through the rest of the gas nearly bal- 
listically. A small fraction of such regions will be on 
nearly radial orbits that impact the innermost parsecs of 
SMBH directly without g oing through a large scale disc. 
INavakshin fc Kingl (J2007D argued that such flows create 
small (pc-scale) warped discs with a fluctuating orien- 
tation. These discs arc able to fe ed the SMBH at much 
higher rates than the flat discs (sec Goodman 2003) with- 
out becoming self-gravitating, effectively due to a much 
faster transfer of gas through the inner parsecs. 

We therefore suggest that the absence of an observed 
Mbh-c relation for pure disc galaxies or pseudobulges 
provides observational support for the chaotic/ballistic 
mode of SMBH feeding. In our interpretation of the data, 
the lack of classical bulges in these galaxies shut down 
the most efficient channel by which SMBH grow - the 
cold chaotic accretion - and therefore led to their SMBHs 
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growing mainly by planar (disk or bar) accretion which 
we argue is inefficient. Having said this, we note that it 
remains to be seen if Mbh ~ & relation for chaotic st reams 
conforms to the relations obtained bv iKind ([2003( 1 . 

6. DEFICIENCIES OF OUR WORK 

We shall now remark on limitations of our work. First 
of all, we only considered the momentum feedback here, 
whereas more complete treatments show that the SMBH 
outflow switches fr om momentum driving to ener gy driv- 
ing at larger radii ([King! 120051: ICiotti et al.ll201fl ). How- 
ever, this is not likely to affect our conclusions signif- 
icantly because it is the momentum-driven part of the 
AGN feedback process that con stitutes the bottl e neck 
for removing the gas to infinity ([King et alj|201lh . Fur- 
thermore, numerical simulations of Zubovas and Nayak- 
shin (2012) that include both the momentum and en- 
ergy feedback forms confirm the points made here on the 
difficulty of expelling a massive dense gas disc for the 
particular case of the Milky Way. 

In common with other analytical papers on the subject 
(e.g.. ISilk fc ReeslfT998t iFabianl fl99fl IKuTbI [2001 [20051: 
iMcLaughlin et al. 2006) we have also used an isothermal 
potential. However the discussion in the end of section 
[2] shows that this has little effect on momentum driven 
outflows studied here. Zubovas & King (2012, in prepa- 
ration) have extended this result to large scale energy- 
driven outflows. 

A more serious worry is that our theoretical model pre- 
dicts Mbh oc c 4 , whereas different authors now find that 
a ranges from somewhat below 4 to as high as 6 (cf. refer- 
ences in the Introduction). One possible explanation for 
this within our momentum-driven feedback model is that 
the steeper slope for samples that combine different types 
of galaxies actually results from a superposition of sev- 
eral Mbh oc <r 4 relations for different galaxies vertically 
offset in mass Zubovas & King (2012, in preparation). 
Future theoretical and observational work is needed to 
resolve this issue. 

Finally, given our SMBH feeding focus, we do not at- 
tempt to relate the gas discs we consider here to expo- 
nential stellar discs observed in real galaxies. This would 
require an additional multi-parameter model for transfor- 
mation of gas into stars and is beyond the scope of our 
paper. 

7. DISCUSSION 

Previous analytical treatments of the M^-a relation 
have modelled the impact of momentum d riven AGN 
feedb ack on a spherical gas distribution (e.g. lKineJl2003l 
12005ft . We have shown analytically that if SMBH is fed 
by a large scale gas disc, there exists a "disc Mbh _ & 
relation" . The latter has th e same shape as the "spheri- 
cal" M bh - o- relation EEi (|2003l 12005( 1 but is offset to 
higher masses by a factor of a few to ten, realistically. 

We emphasise that the precise value of the offset is 
model dependent, but the general trend is undeniable 
even for AGN feedback models not considered her e, such 
as energy-driven models (e.g.. [Silk fe Reeall998[ ). This 
follows sim ply from the geometrical a rgument first em- 
phasised bv lNavakshin fc Power! (12010( 1: a dense flat disc 
and even more so a bar presents a much smaller target 
for AGN feedback to act upon than a quasi-spherical gas 
distribution, and therefore much more massive SMBHs 



could be built if their masses were limited by a feedback 
argument. We may further allude here to the analogy 
from the field of massive star formation, where radiation 
pressure was once thought to limit stellar growt h; in- 
stead, recent simulations bv lKrumholz et all (|2009( 1 show 
that accretion through a disc (which we note is stable to 
self-gravitational instabilities on the scales of the simu- 
lations) is very efficient as radiation is channelled away 
from the disc. If a similar situation applied in galaxy 
formation, SMBHs would grow as long as there is fuel in 
the massive galaxy-wide discs or bars. One could then 
expect SMBH becoming comparable in mass to galaxy 
gas discs themselves, which is clearly obscrvationally not 
the case. 

Reconciliation of this theoretical result with obser- 
vations leads to interesting and significant conclusions 
about th e dominan t mode by which SMBHs grow in 
galaxies. iHul (|2008( 1 showed that pseu dobulge galaxies 
tend to have undermassive SMBHs, and iGrahaml ([2008( 1 
showed that barred galaxies tend to have SMBHs that 
are a factor of ~ 3 to ten smaller than elliptical galaxies 
at the same velocity dispersion. While iKormendv et al.l 
([2011( 1 argue, in variance with the above authors, that 
that SMBHs do not even correlate with velocity disper- 
sions of galaxy discs and pseudobulges, we note that their 
results also show that SMBHs in pseudobulge systems 
are typically undermassive with respect to their classica l 
bulge cousins (e.g., see Fig. 2 in IKormendv et aLll2011( l. 

We suggest that these observations hint that SMBHs in 
most galaxies (1) do not grow by accretion of gas through 
~ Kpc gas discs, most likely because gas is turned into 
stars too efficiently, and (2) grow via a mode directly 
linked to existence of a classical quasi-spherical bulge. In 
this case it is natural that the SMBH in galaxies lacking 
the classical bulge (pseudobulge galaxies) are undermas- 
sive, and it is also natural that SMBHs do not correlate 
with galaxy discs. For the SMBH accretion modes phys- 
ically linked to classic al bulges we suggest the stocha stic 
(|King fc Pringlel2007H or ballistic dHobbs et al.ll201lH ac- 
cretion modes in which the bulge is the source of low and 
fluctuating in direction angular momentum material. 

An additional indirect evidence supporting our conclu- 
sions comes f rom recent sem i -analy tical models of galaxy 
formation by iShankar et al.l ([2012( 1 . These authors show 
that semi-analytical models which turn off SMBH feed- 
ing during pseudobulge growth episodes driven by secular 
instabilities (bars) provide a much better fit to the data 
than models that allow a con current growth of S MBH 
and pseudobulges (see Fig. 1 in lShankar et al.ll2012f ). We 
here thus provide a physical interpretation of this empiri- 
cal result: "fragmentation catastrophe" of the outer cold 
and massive regions of AGN discs makes them too sus- 
ceptible to star formati on, so that not enough gas trick- 
les do wn to the SMBH (|Goodmanll2003HNavakshin et al.l 
l2007tl . 

To avoid confusion, we emphasise that we do not ar- 
gue that small (e.g., sub-parsec) scale accretion discs are 
unimportant for SMBH growth. In the chaotic accretion 
picture, there will always be a non- vanishing amount of 
angular momentum associated with the accretion flow 
and this results in the formation of "small" by galac- 
tic standards discs (R ,$ 0.01 — 0.1 pc, depending on the 
SMBH mass and the accretion rate, cf. iGoodmanl 1200311 
that are resilient to AGN feedback (as we show in the 



Appendix). Rather, our key point is that the origin of 
material that feeds these small discs is distinct from the 
larger scale (with respect to the SMBH's sphere of influ- 
ence) disc or a galactic bar, and the orientation of these 
smaller-scale AGN discs is likely to be chaotic because 
of the manner in which the quasi-spherical inflows that 
feeds these discs fluctuates in time. 

We also note that large scale galaxy- wide gas discs and 
gas bars are obviously important for galaxy formation as 
a whole, likely accounting for a significant fraction of all 



stars formed and also for a radial redistribution of gas 
and stars within the galaxy, but we see no clear evidence 
that these features are able to channel their gas all the 
way to the SMBHs. 
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APPENDIX. FEEDBACK RADIUS. 

We have argued above that AGN feedback can affect a large scale gas disc expelling the gas outward just as it can 
for a spherical distribution of gas. We now show that there exists a critical radius within which this argument no 
longer holds: the disc there is resilient to AGN feedback. 

Inside the SMBH's sphere of influence, marked by the radius i?/> = G M^/a 2 , the gravitational potential is domi- 
nated by the SMBH and so we can replace M(R) with Afbh- For a SMBH on the Mhh~<J relation, the SMBH's radius 
of influence is given by 

R h « 10 pc M 8 a 20 2 . (11) 

As the result, the critical momentum flux striking the edge of the disc must be 

R -f* * -V M <- ■ < 12 > 

Using A/disc ~ Mbh(-ff/-R) for a marginally self-gravitating disc again, we arrive at the requirement that 

k£ bh - k^§ < 4 • (13) 

nR z 

This can be satisfied at radii greater than 

R> (KM bh /47r) 1/2 «25pcM 8 1/2 , (14) 

where M§ is the SMBH mass in units of 10 8 M Q . Recalling that we assumed that R < Rh, we conclude that inside 
the "feedback radius" _Rfb given by 



-Rfb 



i? ?l ,(KM bh /47r) 1 / 2 , (15) 



momentum feedback by the SMBH is unable to affect a self-gravitating disc. 

Physically the significance of this feedback radius can be appreciated by considering a major gas mass deposition 
event - say, the infall of a very massive Giant Molecular Cloud or a merger with a gas-rich galaxy that dumps a 
significant amount of gas into the inner galaxy. This gas will settle into a disc because it is very likely to have a non 
zero net specific angular momentum. At R <C Rr,, the disc cannot be expelled by the AGN, whereas at R <C R^, it 
may as shown in the main body of the paper. 

Somewhat counter-intuitively, gas closest to the SMBH is actually hardest to expel to infinity (although this becomes 
clearer if one considers that SMBH's gravity is enormous in the inner ~pc). Deposition of a small cloud at R <?C i?fb 
may be much more promising for AGN feeding than having a vast amount of gas at R ^> R^ as the latter is more 
easily affected by AGN feedback and star formation. 



